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We investigate higher twist contributions to the transverse momentum broadening of Drell Yan 
pairs in proton nucleus collisions. We revisit the contribution of matrix elements of twist-4 and 
generalize this to matrix elements of arbitrary twist. An estimate of the maximal nuclear broadening 
effect is derived. A model for nuclear enhanced matrix elements of arbitrary twist allows us to give 
the result of a resummation of all twists in closed form. Subleading corrections to the maximal 
broadening are discussed qualitatively. 
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I. INTRODUCTION 

In the era of RHIC and LHC the understanding of 
hard processes in heavy ion collisions plays an increas- 
ingly important role. It has become clear that the onset 
of the regime of perturbative quantum chromodynamics 
(pQCD) with increasing momentum transfer in reactions 
involving large nuclei can be very different from that in 
collisions of single hadrons. Previous studies have shown 
that the leading twist approximation |Q] in pQCD, even 
at intermediate momentum transfers, can be spoiled by 
higher components of the twist expansion. 

The reason for that is that matrix elements of higher 
twist can be increasingly sensitive to the size of the nu- 
cleus. With naive power counting, twist-2 parton distri- 
butions scale with the mass number A of the nucleus, 
while matrix elements of twist-(2m + 2), m G N, can 
scale with m additional powers of A^/"^. These matrix 
elements are then called nuclear enhanced. They corre- 
spond to processes with multiple scattering. For large 
nuclei these additional factors of A^^^ can compensate 
the inherent power suppression of higher twist. 

Let us explore the physical picture Q behind this for 
the example of p -I- ^ collisions. In the leading twist ap- 
proximation only one parton a from the nucleus and one 
parton b from the proton participate in the hard scat- 
tering a + b. They are described by parton distributions 
respectively. A nuclear enhanced twist-4 matrix element 
for the nucleus corresponds to two partons a and a' en- 
tering the hard scattering, They enable us to describe 
the process where the parton from the proton scatters 
off these two partons: aa' + b. The additional factor 
of A^/^ arises when the two partons come from different 
nucleons inside the nucleus. This mechanism was first 
pointed out by Luo, Qiu and Sterman some years ago 

^, P based on generalized factorization theorems 

In the following years various calculations on the level 
of nuclear enhanced twist-4 (double scattering) became 
available for the Drell Yan process p+A in proton 

nucleus collisions [0, ^, ^- However going beyond the 
level of twist-4 remained a technically difficult point. It 
is known that for the Drell Yan process in symmetric 
colliding systems the QCD factorization fails at the level 
of twist-6 ]lO| . However iTip + A we will formally work in 
the limit of very large A and therefore stay in the leading 



twist approximation for the proton. This allows us to go 
to arbitrary nuclear enhanced twist for the nucleus 
Nevertheless a systematic calculation of contributions of 
arbitrary twist was never done in this framework. 

In this publication we study the transverse spectrum 
of Drell Yan pairs and its nuclear broadening. This was 
already done on the level of twist-4 We will repeat 
this calculation both in light cone and covariant gauge 
for the strong interaction. In particular we comment on 
the question of electromagnetic gauge invariance. After 
that we generalize the result by calculating the contribu- 
tions of operators of arbitrary twist. This is feasible for 
contributions which lead to maximal broadening, there- 
fore establishing an upper bound on the nuclear effects on 
the transverse momentum spectrum. Besides this there 
are subleading contributions which we will also discuss 
shortly. We then present some arguments how the higher 
twist matrix elements can be approximated by model de- 
scriptions in such a way that we are able to give a closed 
expression for the sum of all higher twist contributions. 



II. NUCLEAR ENHANCED DOUBLE 
SCATTERING 

In this section we will go through the calculation of the 
twist-4 contribution to the nuclear enhanced Drell Yan 
process. This computation was already done in covariant 
gauge by X. Guo some time ago We confirm her re- 
sults. We also present an explicit calculation in light cone 
gauge, which bears some difficulties for higher twist. We 
also comment on the problem of electromagnetic gauge 
invariance in this calculation. 



A. Light cone gauge 

The Drell Yan cross section, integrated over the an- 
gular distribution and the rapidity of the lepton pair, is 



da 



E 



247r2S'Q2 



(1) 



where Q, q± and y are mass, transverse momentum and 
rapidity of the virtual photon and S = {Pi+P2)'^- Pi and 



2 



P2 denote the four momenta of the nucleus (per nucleon) 
and the proton in the center of mass frame respectively. 
The sum runs over all quark flavors a and Ca is the charge 
of a given quark flavour relative to the electron charge. 
The hadronic tensor is defined as the matrix element of 
two electromagnetic currents 

W^. = J d^ye^^-y (P1P2I jnO)r(y) \P1P2) (2) 

in the colliding system. After integrating the angular 
distribution the tensor P^'' reduces to the polarization 
tensor of the virtual photon 

P"" = -g'" + ^ (3) 

where q is the four momentum of the virtual photon 
Electromagnetic gauge invariance of course requires that 



q'^W^,, = = W^,q\ 

We want to calculate the nuclear enhanced twist-4 con- 
tribution to the hadronic tensor. Thence we have to take 
into account the diagrams shown in Fig. |^. We denote the 
momenta of the quarks with ri and r2 and the momenta 
of the gluons with K and K' . We expand the momenta 
of the parton lines between the hard and the soft part 
into a longitudinal contribution on the light cone and a 
transverse contribution 

< = eiPr + r^ r^2=^2P^ + ri (4) 

= a;Pf + K^t" K'f' = x'Pf + K'^" (5) 

introducing parton momentum fractions x, x' , ^1 and ^2- 
The hadronic tensor for the symmetric diagram (a) is 



Anas 

1 



dy 



d^d^ 
6 6 



P+dxP^dx'd'^rrtd'^K^d'^K'^fq{^2){S2TT^)T'-''^Hj^J^ 



"•"'P" (rs + )2 + ie 



(r2 + K' 



dy, d^yj dy^ d^y^ dy^ d^^yi ^-^i^P+v7 .-^xP^V^ .^x'P+v7 ^-^r^■V^ ^~^K^-V^ ^^K'^■v^ 

(2^)3 (27r)3 (27r)3 



1 



{Pi\mi^q{yi,vi)A'{v^.yi)A''{vl.yi) lA) (6) 



for the case that the quark comes from the nucleus and 
the antiquark comes from the nucleon. Here 



x5[ql + {ri +K^f]S 



1, a + 

y — In — T- 
2 6 



(7) 



are the momentum conserving 5-functions which connect 
the photon momentum to the internal momenta. It is 
understood that the (5-functions have to be integrated 
in order to obtain physical results. We have chosen to 
integrate the rapidity dependence of the cross section 
which is trivially done by the last (5-function. The first 
(5- function will be eliminated by the integral over ^2- At 
a later stage we will also integrate over q\ but since we 
will do that in different ways it is a matter of convenience 
to keep the unphysical J-function in the mean time, fq in 
above formula is the antiquark distribution in the proton. 

The hard part TL of the cross section is given by the 
propagators in Eq. (||), the phase factor F^"^ and the 
trace 



/ 2 respectively. In order to project out the correct twist 
contribution there is an important step. We have to per- 
form a coUinear expansion of all momenta entering the 
hard part which is equivalent to a Taylor expansion in all 
transverse momenta. In light cone gauge for the gluon 
field we can stick to the 0th order of the expansion, i.e. 
we simply set all transverse momenta equal to zero. In 
Eq. (^) this was already done for the quark lines from 
the nucleon since we only wanted to generate the ordi- 
nary twist-2 parton distribution fq{£,2)- For parton lines 
from the nucleus we keep the transverse momenta for 
the moment. The reason is that this coUinear limit is 
not straightforward in light cone gauge. In fact, it will 
turn out that some contributions which are falsely taken 
to vanish in the coUinear limit in light cone gauge give 
finite results because the zeros in the hard part are can- 
celed by poles elsewhere. 

The two remaining diagrams (b) and (c) in Fig. ^ have 
the same decomposition of the hadronic tensor but with 
denominators of the propagators replaced by 



H\^lp. - ^Tr[7.(-/2-^')7./27p(-/2-^)7M/i]- (8) 

The quark lines connecting the hard part to the nucleus 
and the nucleon respectively are contracted with f i and 



1 



1 



(r2 +K' - KY - ie (rz + K' f - ie 
1 1 



{r2+K- K'Y + ie (ra + + ie 



(9) 
(10) 
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FIG. 1: The diagrams contributing to Drell Yan on the level of twist-4. The symmetric diagram (a) is double scattering, the 
asymmetric diagrams (b) and (c) correspond to an interference of single and triple scattering. 



for diagrams (b), (c) respectively, 5- functions 



xS[ql+rn6 



2 6, 



(11) 



X 5 



y in 

2 6 , 



and traces 



(12) 

(13) 
(14) 



Next we want to convert the gluon gauge fields to field 
strengths. In light cone gauge = the transverse 



components of the A fields give the dominant contribu- 
tion and we can use 



P+APP+A" = —{d+AP)-^[d+A'') 
ix —ix' 

«^F+-F^+(-gf/2) (15) 
xx' 

By the introduction of the field strength tensors we make 
the poles in the parton momentum fractions x and x' 
explicit that were hidden in the gluon gauge fields. It 
will turn out that these poles are artificial and will be 
canceled by zeros which are hidden in the hard part. This 
is a special feature of the calculation in light cone gauge. 

We carry out the coUinear expansion as the lowest or- 
der Taylor expansion of the hard part in the transverse 
momenta , and K'^ 

n{±) = H(±)|^^„+^ ± _^H(±)|^^„+o(±). (16) 

Here _L shortly stands for all the transverse momenta. 
Now we are able to carry out all transverse integrals and 
to eliminate all dependences on transverse coordinates in 
Eq. (g) in a trivial way. For diagram (a) this gives us 



jdyj j dxdx'm^)^ [H{ri,K\K'^)]^^^{-g'r/2) 

J ^^^y4-^-^«^'^^"^^"e-"^^"^3-e-'^."^.-i (Pi|,-(0)7+g(yr)^^+"(% )^.+(y^ \Pi) 



(17) 



Due to the subtle interplay of poles and zeros in light 
cone gauge we evaluate the hard part with transverse mo- 
menta however keeping in mind to take the limit _L— > at 
the end. There are some caveats. E.g. if we want to test 
electromagnetic gauge invariance it is clear that we must 



not take q^^ limj^^o W^iy but we must take limj^_+o q^'W^i, 
since the hadronic tensor contains (5-functions which give 
a finite transverse component as long as the internal 
transverse momenta are still finite. 

The denominators of the propagators can be trans- 
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formed into poles in x and x' via 
1 1 



(r2 + K)^ + it (r2 + K'Y - ie 



(18) 



We define the prefactor and the tensor structure from the 
field strength into the trace and finally evaluate 

Hi'i) = {^2S)-'Hl^l, (-.9f /2) 
6 



6-5 

1 



2xx'Pi ^Pi , + xPi + x'P^ 

1 



■ K' 



-iPlf.P2.+P2^Plu) 



(19) 



At this stage we see that the naive limit ±— > carries 
most of the terms in above equation to zero. However 
if we look at it more carefully through the factor 1/xx' 
from the gluon fields we will encounter a 0/0 situation 
which we have to resolve in order to obtain finite results. 

For diagrams (b), (c) we get similar poles and the 
traces reduce to 

- iC2S)-'Hl%^ (-gf /2) 



6^ 



Pi,iK'^-K^), 



S,2{Pl fj,P2 y + P2 iJ.Pl p) — -^^2Sgp,i, 



(20) 



H'i} = {C2S)-'Hl^l^ (-.9^/2) 



^2S 



P,^{K^-K'^), 



UPl^P2, 



P2 tJ.Pi u) - -^^2Sgfji, 



(21) 



We carry out the integrals in x and x' by taking the 
residues given by the poles from the propagators. By the 
usual trick we have to take into account the exponential 
factors in the matrix element to argue in which way to 
close the integration contour and thereby we introduce 
an ordering of the light cone coordinates by 



0(c) = 



0(yr-%)0(-2/4"), 
0(yr -2/^)0(2/^ - % ) 



(22) 
(23) 
(24) 



for diagrams (a), (b) and (c) respectively. 

At this stage we see that the limits r^,r'i — *■ 0, are 
trivial in the sense that there are no poles or zeros ap- 
pearing by taking this limit. We can therefore set these 
transverse momenta to zero. This implies K'-^ = K-^ 
and x' = X = —K^ /£,2S. Now we add up the contribu- 
tions of all three graphs. Note that the pole integrations 
introduce minus signs for graphs (2) and (3): 



Q(a)^(a)p(i) _ Q(b)H(j>)T(^) - e(^)i/(t^r(3) 



X^l 

6^ 



2xPi^Pi, + {Pi^K^ + Pi.K^) 



0(i)p(a) 



^ [Isg,. Pi,P2. - P1.P2,] ef"^) [rt'^) - rf)]] 



^ilsg,. Pi,P2. Pi.P2M]rW[e('') - e('') - e(^)]. (25) 



Let us now discuss electromagnetic gauge invariance. 
After taking the poles all quark lines attaching to an 
arbitrary quark photon vertex in either diagram are on 
the mass shell. Therefore electromagnetic gauge invari- 
ance at this stage is ensured for each diagram separately 
by simple equations of motion. Let us denote the three 
terms on the three different lines on the right hand side 
of Eq. (^5|) by A, B and C starting from above. It 
is easy to check that C is separately gauge invariant 
and so is the sum of A and B. B and A alone are not 



gauge invariant as the contraction of A with q^^ gives 
1/2 x^iK'^Q^°-'>T^°-^ which is then canceled by the same 
term with different sign originating from contraction of 
with B. As already discussed the 5-functions deter- 
mine = (^1 + x)P^ + ^2P2 + (for r^'')) and 
q>' ^ ^iP^ + £.2P!^ (for r('')). 

We will now drop the term C because the term O^"-* — 
requires the partons from the nucleus to be 
bound together and leads to a twist-4 contribution which 
is not nuclear enhanced pi 0] . The remaining term A+B 
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is only proportional to 0^°) and is manifestly gauge in- variant. We absorb 8^°-' into the matrix element 



2n 27r 



iy^-vI)Qia) (p,\q^O)^+q(^y-)F+-^y-)F^+{y-) \P,) , 



(26) 



We can write the hadronic tensor after summing over 
all diagrams and integrating over and y as 

W,. = ^^^^ f^h{^2) lim T,,{i,,x)H,, (27) 

where B = Q'^/S, ^ = QV(6'5') and 



H^^ — 



[S{ql-kl)-S{ql)] 



(28) 



We have introduced A:^ = -K^ > 0. 

We note that this expression still has poles due to our 
choice of gauge for the gluon fields, but as checked above 
electromagnetic gauge invariance is manifest. It is once 
more interesting to note that the term in above equation 
which is proportional to the transverse gluon momen- 
tum is needed to ensure gauge invariance for the entire 
hadronic tensor. This now justifies the careful treatment 
of the coUinear limit. 

We are ready to contract the hadronic tensor with the 
photon polarization tensor. Since q^H^j^i, is zero this is 
simply a contraction with —g'^'^. The result is simple and 
finally allows to take the limit — > unambiguously, 
leading to 

i?^.P^''= lim ^[Siql-kl)-Siql)]=-5'iql). (29) 
The final result for the cross section is 



da I 



dQ^dqj 



SNcQ^S 



/'f^r,,(6)/,(6) 



Nr. 



(30) 



with the known soft hard matrix clement Tqg(^i) — 
Tqgi^iiX = 0). A second term proportional to Tqgfq 
has to be added but is not explicitly shown. This in- 
deed coincides with the result in We introduced a 
subscript "1" with the cross section to indicate that this 
is the contribution with one pair of gluon fields in the 
matrix element. 



B. Covariant gauge 

Let us briefly go through the corresponding calcula- 
tion in covariant gauge, since this will be the technique 
we use later for arbitrary twist. We can again start from 
Eq. (^ and proceed along the same way as before, keep- 
ing transverse momenta for the gluon fields. However in 
covariant gauge we use 



P+APP+A" « PPA+P^A+. 
We therefore define 



(31) 



(32) 



where HliO,p<T are the traces for the graphs i = a, b and 
c. In covariant gauge we have to carry out the coUincar 
expansion to the second order in transverse momentum, 
applying the operator 



2! ^ dK^dK^ 



(33) 



to the hard part. The explicit factors of transverse mo- 
mentum are used to obtain field strength tensors in the 
matrix element via 

KiA+{y:,)KiA+{y,) 

^F+-iy,)F^ + iy,){-gt/2). (34) 

After performing the pole integrations in x and x' and 
integrating over and y the hadronic tensor for diagram 
(a) is 



-^Iq{£.2)-T- 



9L 



2! dKj^dKi 



Tqg{^ux)H\^j5{ql+K' 



(35) 



where now x = x' = k\/^2S and — Q^/£,2S. Since 

-g'^'Hl:) = 66^, (36) 

the derivatives can only act on the transverse (5-function 
and the phases in the matrix element. 

For diagrams (b) and (c) Eq. (^5|) holds with an addi- 
tional "— " sign and modified G-functions 0^^^ and 0^^^. 
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Furthermore the transverse (5-function is just and 

-g'^'H'i'} = = a^S'. Therefore from the 

asymmetric diagrams (b) and (c) we can only get nuclear 
enhanced contributions with derivatives on the phases in 
the matrix element. When we sum over all diagrams the 
terms with derivatives on phases are again proportional 
to 



Q(a _ Q(b) _ q(c) 







(37) 



and therefore show no nuclear enhancement — confer 
the discussion below Eq. We conclude that the 

only nuclear enhanced contribution can arise from the 
transverse (5-function of the symmetric diagram. We have 



5a. 
2 



This implies 



1 



da I 



2! dKj^dK^ 



47rc^a^ 



-5'{ql). (38) 



dQ'^dq\ SN^Q^S 



/'^T,,(a)/,(6) 

Jb ^2 



Nr. 



5'{q 



(39) 



which coincides with the result in light cone gauge. 

Let us shortly discuss the result of the twist-4 calcu- 
lation, ^-functions and their derivatives, strictly speak- 
ing, are only defined when they are integrated. This ap- 
plies also to our result. We have to integrate over the 
transverse momentum in order to obtain an observable 
quantity. Let us define moments Q^™ of the transverse 
momentum spectrum by 



2 „2m 



dqi q_ 



da 



dQ-^dq^. 



(40) 



We recall the leading twist result for the cross section in 
leading order in as 



dao 



Altera 



' ^Mil)M^2)Siql). (41) 



dQ^dqi SNcQ^SJb 6 

A simple integration over gives zero which implies 
that there is no nuclear enhanced twist-4 contribution 
from these graphs to the spectrum integrated over the 
transverse momentum 



= 0. 



(42) 



The subscript "1" on the left hand side refers to the con- 
tribution from CTi. This is in disagreement with the re- 
sults of However, when we look at the second mo- 
ment of the transverse momentum spectrum, we get a 
nuclear enhanced twist-4 contribution 



2 47r c^Q!^ 



'd^2 



(43) 



From leading twist we obtain 

rn" ^ - ^^^'9"' f' d^2 



/?(6)/g(6 



{Ql)o = 



(44) 
(45) 



for the first two moments. All higher moments vanish 
both for twist-2 and twist-4. 

The soft-hard matrix element is usually approximated 

by 



r,s(a)«A2^4/V,(6) 



(46) 



where fq is the nuclear parton distribution for a quark 
normalized to one nucleon and = 0.01 GeV^ @, |, [|. 
This leads to an additional broadening of the transverse 
momentum spectrum in p+A collisions compared to p+p 
collisions 



A(gi) 



[Qj_]pA 



pA 



[Q±]pp 



(47) 



pp 



To twist-4 accuracy we have Qf' = (Qi")o + (Qi")i- 
Using the model for the soft-hard matrix element we ob- 
tain Q 



(48) 



The transverse momentum broadening with the charac- 
teristic A^^^ scaling was observed by the Fermilab E772 
collaboration IH. 



III. BEYOND DOUBLE SCATTERING 
A. Preliminaries 

The result of the last section raises the question 
whether terms of even higher twist contribute to higher 
moments of the transverse momentum spectrum. This 
is indeed the case and will eventually allow us to resum 
the contributions of different twist. To start we want to 
pin down what we are looking for. We are interested in 
contributions of Twist-(2n + 2), n e N, with maximal 
nuclear enhancement A^/'^. Therefore we take into ac- 
count operators with n pairs of soft gluons, each single 
pair a color singlet with no flow of momentum between 
these pairs. These are necessary conditions to achieve 
maximal nuclear enhancement. The diagram in Fig. |^ 
is an example with 2n gluons. Furthermore we use the 
argument, that in a large expansion the leading con- 
tribution is coming from planar diagrams [p^ . In order 
to handle the large number of possibilities to group 2n 
gluons together to pairs we stick to the leading order in 
1/A^c and only take into account planar diagrams, i.e. 
with no gluon lines intersecting each other. See Fig. ^ 
for an example 

We have seen that the case of twist-4 was simpler to 
handle in covariant gauge. Therefore we will use this 
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FIG. 2: Scattering off 2n gluons. In this example k gluons 
(k < 2n) are on the left side and 2n ~ k gluons on the right 
side of the cut. The gluons are ordered to pairs with momenta 
Ki and K'^ for one pair. Big dots indicate soft gluon poles. 



gauge here. This requires us to take a derivative of the 
hard part of order 2n with respect to the transverse mo- 



menta. We denote the gluon momenta with Ki, K2, ■ ■ ■ 
and K[, K!^, ■ ■ ■ , where always Ki and K- belong to one 
pair. We introduce parton momentum fractions by 



K' ^x'Pi+K 



(49) 



The differential operator we have to apply is of the form 



(2n)l 



dK^^^dK'^"^ 



(50) 



In order to obtain the matrix element with 2n field 
strength tensors in covariant gauge, we have to apply 
derivatives to all 2n momenta once. To keep it simple 
we investigate the contribution with a maximal number 
of derivatives on the transverse (5-function. This will give 
us a maximal broadening effect. 

E.g. the diagram in Fig. || provides a term S{q'j_ + {K^ + 
. . .+K^)'^). When we keep our requirement in mind that 
there should be no flow of momentum between the pairs 
of gluons we have K^ = K[^ . Then application of the 
operator (pG) to this transverse 5-function gives 



J 



J2n 



{2n)\Y[^dK^^^dK^^^ 



5{qi + {Kt + --- + KiY)^5.. 



2" 



.<5(")(,i)n/*^' 



(51) 



r 



where denotes a nth derivative of a (5-function. 

In the complete Taylor expansion of order 2n there are 
(2n)!/2" derivative terms of the kind of Eq. ( [50| ) which 
have to be summed. This exactly cancels the combina- 
torial factor in above equation. We note that we only 
get a non-vanishing contribution if fc = n, i.e. if the di- 
agram is symmetric. It is straightforward to show that 
for fixed n the symmetric diagram in Fig. |^ is the only 
planar diagram with 2n gluons which contributes with 
the maximal number of transverse derivatives. Therefore 



with the restrictions we have chosen we have only one 
diagram to calculate in covariant gauge for each n. 



B. Calculation of arbitrary twist 



The hadronic tensor for the diagram with n gluon pairs 
can be written as 



J 



W" = (4™,)"C„ 



d^d^ 



Yl dxMd^K^d^Kl^ /,-(6)i?;,(32^3)r"p- 



+„- T-r / dy- dy', _,-c.d+„- „„'d+„'-1 



2ti 2ti 



e~^i^P^y^ e^<P^y'^' Up^\mi^q{y,)\{A+{yi)A+{y'r) \P,) . (52) 



We start with the evaluation of the color factor for every For the hard part we get the simple result 
diagram which is 
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{-gn\{Pi'P^' =U2S. (54) 



This can be proved by repeatedly applying the identity 

fi(...)fi=2Pi.(...)f i=6^fi (55) 

where (...) stand for one of the parentheses in the equa- 
tion above. 

The 2n propagators provide poles 



7?" 



1 



1 



Xi 



52 s 



+ ie xi + X2 + 



«2S 



■ ze 



(56) 



and we also have 5-functions 



x^('/i + (E^'0') 



X (5 ( y - 2 In — ^^"^ 



(57) 



After integrating the momentum fractions Xi and x[ by 
taking the poles and integrating over y and we apply 
the operators from Eq. (|5^) to the transverse ^-functions 
and use the explicit factors of if^^ in front of the operator 
together with Eq. (53). We obtain a hadronic tensor 



327r^ 



X / --/,(6)T,,.(a)'5^"Hgi) (58) 

where = /S,2S and the matrix element with n pairs 
of gluon operators is defined as 



2tt 



n 



27r 



X 2(^il9(0)7+9(yo") 



xn^^''(2^r)^^ + (2/rM%ni^i)- (59) 

The 0-functions from the pole integrations reflect causal- 
ity and give an ordering of the vertices along the quark 



lines on the light cone 

6" - Q{y^ - y-)e{y- - y-_,) ■ ■ ■ Q{y^ - y^) 
X Q[-y'-)Q{y'- ~ y'-^,) ■ ■ ■ e{y'f ~ y'f). (60) 
The cross section can finally be written as 



dQ^dql SNcSQ^ 



47r^Q!, d 



dql 



I f^T,,.(a)/,(6), (61) 



including the known cases for n = 0, 1. Again we can 
obtain physical results only if we integrate over the trans- 
verse momentum and resolve the 5-functions. From ma- 
trix elements with n gluon pairs we get contributions to 
the 2mth moment of the transverse momentum spectrum 



''SNcSQ^ 



'dC2 

< 6 



T,gn(a)/,(6)- (62) 



When applying the derivatives in the differential op- 
erator ( |50| ) to the hard part we have chosen to take all 
derivatives on the transverse (5-function. We know that in 
the twist-4 case (n = 1) this was the only way to obtain 
nuclear enhancement. This is no longer guaranteed for 
n > I. The derivatives can also act on the phases in the 
matrix element. Moreover more graphs and orderings of 
gluon pairs can contribute. This has to be taken into ac- 
count for a more quantitative analysis. In general, from 
a diagram with n gluon pairs we expect a cross section 
which can be written as a sum 



da 



dQ^dqj 



E 

fc=0 



Vk 



Q2(n-fe) 



jW(,i). 



(63) 



Lower derivatives in are compensated by explicit fac- 
tors of The Vk are coefficients which are unknown 
at present except for the case k = n where they can be 
taken from Eq. (|6l|). 

Therefore from a diagram with n gluon pairs we would 
get terms contributing to all momenta of the spectrum 
lower than n. We were already arguing that our result, 
where we only keep the terms k = n, will give us a max- 
imal broadening effect. We will discuss later how this 
shows up after summation of all twists and how the sub- 
leading terms with less derivatives will alter the result. 
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FIG. 3: The symmetric diagram with n gluons on each side, 
where opposite gluons constitute pairs of color singlets, is the 
only planar diagram which contributes to maximal nuclear en- 
hancement. Dotted lines show the planar flow of color through 
the matrix element. 



C. Resummation 

In order to sum contributions of arbitrary twist, one 
has to establish relations between the matrix elements 
Tqgn of different twist. Our approach will be to develop 
a simple model for these matrix elements. The usual 
way to treat the twist-4 soft hard matrix element is to 
factorize the soft gluons from the quarks. For that one 
inserts complete sets of nuclear states 



(64) 



between the pairs of parton operators. A'' is a normal- 
ization factor. Of is a symbolic index for all the excited 
states of the nuclear system and p is the momentum nor- 
malized to one nucleon. Later we approximate sum and 
integral by the contribution of the lowest state with the 
given momentum We assume that this is a good 

approximation when there is no momentum transfer be- 
tween the parton pairs. 

Furthermore we introduce relative and absolute coor- 
dinates for the gluon pairs 



(65) 
(66) 



Since gluon pairs are bound to color singlets we can as- 
sume that ^ Ra for large A, where Ra is the nuclear 
radius. This enables us to use 

6" « e(-Z„)e(Z„ - Z„_i) • • ■ 6(^2 - Zi). (67) 

Corrections to that will lead to contributions which are 
suppressed by powers of A^/'^ compared to the full nuclear 
enhancement. 



We introduce bilocal gluon correlators 

dz. 



G{Z, 



2nN 



(Pi|F+'^(Z 



z,)F^ + {Z, + z,)\P,). 

(68) 



We assume now that the value of this correlator does not 
depend strongly on the absolute position of the gluon pair 
inside the nucleus, so that G{Zi) « Gq. Our assumptions 
lead to 



r,,„(0 = ^/,(OGJ? / \[dz, 

X 0(-Z„)0(Z„ - Z„_i) • • • 8(^2 - Zi) 
= Am)-,\{(Gj dZ, 



(69) 



The term in parentheses scales with A^/^ and we set it 
equal to X'^A^/^. was already introduced for modeling 
the twist-4 matrix element. Finally we have 



(70) 



As expected the twist- (2n+2) matrix element scales with 
n additional powers of A'^/'^ . 

The assumptions we have used can be summarized by 
the statement that all QCD dynamics is reduced to the 
QCD substructure of the single nucleons in the nucleus, 
whereas the nuclear effects are taken to be purely ge- 
ometrical in nature. This might be quite reasonable to 
estimate leading effects. We remember that in the case of 
parton distributions the idea of a superposition of parton 
distributions of the individual nucleons is also leading to 
a reasonable first guess for the nuclear parton distribu- 
tions. Dynamical nuclear effects then enter as corrections 
which might be as large as 20 or 30%, showing up as phe- 
nomena like shadowing and the EMC effect . 

When we plug our model into Eq. (61), we get 



dun 



.2^,2 



dQ^dqj 



'SNcSQ^Jb 6 
1 



mi)m2 



Nr. 



dq^ 



S{ql). (71) 



Summing up all twists leaves us with a shift operator, 
acting on the 5-function 



da 



dQ^dq^^ 
= exp 



E 



da„ 



dQ^dq^ 



Nr. 



l!l!f!f./|l/3\2_^L 
dql 

d(To ( 2 



dQ'^dql 



Nr 



A'/^X"^ (72) 



At first glance this is a surprising result. Let us note 
that the summation of an infinite series of derivatives 



10 



of ^-functions can give a well-defined smooth function. 
That our particular result is again of S shape is in some 
sense accidental and due to our simple model. 

However we remember that we already expected a re- 
sult which gives us a maximal nuclear broadening. This 
was because we dropped all terms with less than n trans- 
verse derivatives in the general expression in Eq. (|63|). 
The result of a shifted (5-function, where there is no 
contribution with m — left, is the outcome of cut- 
ting the sum in (|63|). We expect contributions with 
less derivatives on the transverse (5-function to alter the 
shape through explicit dependence which is absent 

in our result. This should lead to contributions which 
are less shifted and therefore, when added, filling the 
"valley" between q± = and the maximal shift value of 
ql = ATT^asA^^^XyNc. 

A more rigorous argument would be that for a given n, 
because of the lower number of derivatives in the terms 
dropped in Eq. (|6^), these add to lower moments, there- 
fore favouring lower values of in the function recon- 
structed from these moments. A more quantitative anal- 
ysis is needed here. We shall report on that in a forth- 
coming publication. However our result allows to give the 
order of the smearing of transverse momentum through 
nuclear effects. It is of the order of 

ql = An^asA^^^X^/Nc ~ (400 MeV)^ (73) 

for large nuclei and for A'^ = 0.01 GeV^. 

We have to keep in mind that we are dealing only with 
additional nuclear effects that add to the smooth trans- 
verse momentum spectrum in p + p collisions and which 
can be described theoretically by higher order calcula- 
tions and resummation of radiative corrections. A re- 
summation of nuclear higher twist is an additional effect, 
though not completely independent from a resummation 
of radiative corrections. The scale of these nuclear effects 
is determined by q± . This was already claimed to be the 
scale some time ago based on pure twist-4 calculations [|| , 
but is now confirmed to hold also if one takes into account 



matrix elements of arbitrary twist. q± has also been con- 
nected to the work of Baier, Dokshitzer, MuUer, Peigne 
and Schiff on nuclear momentum broadening [ p7[ . In this 
work also multiple scatterings are effectively summed up. 



IV. CONCLUSIONS 

We have calculated the transverse momentum broad- 
ening of Drell Yan pairs in p + A collisions by taking 
into account nuclear enhanced higher twist effects. We 
presented explicit computations of the twist-4 contribu- 
tion in light cone gauge and covariant gauge, confirming 
earlier results in covariant gauge. 

Furthermore we obtained results for the maximal 
broadening effect of leading order diagrams with arbi- 
trary twist. By relating the corresponding higher twist 
matrix elements through a particular model we were able 
for the first time to give a closed form for the sum over all 
twist contributions. We found that the effect of nuclear 
broadening is of the order of q± ~ 400 MeV for large 
nuclei. 

The exact shape of the broadening effect can be ob- 
tained by the study of terms which are subleading in the 
number of derivatives in transverse momentum. This has 
to be investigated in the future. 
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